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Abstract: In the large extra dimensional ADD scenario, Z bosons undergo a one-loop 
decay into a photon and Kaluza-Klein towers of gravitons/gravi-scalars. We calculate 
such a decay width, extending previous arguments about the general form of the four- 
dimensional on-shell amplitude. The amplitudes calculated are relevant to processes in 
other extra dimensional models where the Standard Model fields are confined to a 4-brane. 
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1. Introduction 

Nieves and Pal have considered the decay of a Z boson to a photon and a standard 
four-dimensional graviton. In the ADD scenario |^]^, the graviton may be viewed from a 
four-dimensional perspective as gaining a "tower" of massive Kaluza-Klein (KK) excita- 
tions m, ^ (this tower takes a relatively simple form if we assume that the extra dimensions 
are toroidally compactified with the copies of having a common radius R/2tt, with R 
the common circumference). Although the decays of a Z boson involving real production 
of a KK graviton excitation with a photon will be suppressed by a gravitational coupling, 
the existence of a "tower" of particles to which the Z boson can decay may counteract this 
suppression. Because the KK excitations only couple with gravitational strength to Stan- 
dard Model particles, they will almost certainly pass through a detector (their detection is 
a next-to- leading-order process in the gravitational coupling). The ADD scenario therefore 
predicts that we should see an increase in the decay width of the Z boson to a photon and 
missing energy relative to the Standard Model prediction |^ |^, |ll| . 

In this paper, we calculate, to leading order, each of the decay widths of the Z boson to 
a photon and a KK graviton excitation. There are two relevant towers of KK excitations: 
a spin-2 tower and a spin-0 tower. We combine the calculated widths to obtain an overall 
decay width for the decay of a Z boson to a photon and some KK graviton excitation. 
This calculated width will allow the determination of bounds on the size of large extra 



dimensions when combined with experimental data on Z decay |12|. (The current upper 
limit on the branching ratio for Z — > 7 + X, with X some beyond-Standard Model invisible 
particle or particles, is O(10~^) ||l^. With a "Giga-Z" collider setup, this could potentially 
be reduced to around O(10~^).) The amplitudes calculated are also relevant to processes in 
other extra dimensional models where the Standard Model fields are confined to a 4-brane 
(e.g. the Randall-Sundrum 1 (RSI) model 0]). 

Current experimental limits on the size of ADD extra dimensions from processes other 
than Z boson decay come from inverse square law experiments |]l^, 16|, from consideration 



of the channel e"*" + e 7 -|- graviton at the LEP experiments [17|, and from consideration 



^The idea of large extra dimensions was also considered ^ ^ ^ prior to the work of Arkani-Hamed, 
Dimopoulos and Dvali. 
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of the channel p+p — > jet+graviton at the Tevatron experiments [18, 19]. In the near future, 
the most likely improvement in these experimental bounds should come with the publication 
of results from the D0 experiment using Run II of the Tevatron, and combination of those 
results with the already published CDF Run II bounds Further into the future, we 
can expect investigation of large extra dimensions at the Large Hadron Collider. 

For the decays considered in this paper, the leading order process is a one-loop process, 
because the tree-level vertex for Z boson decay into a photon and a graviton is absent. This 
means that there is a higher order prefactor in coupling constants of the decay width when 
this decay mode is compared to other, tree-level, graviton production modes considered 
previously 0, |2^. However, there is a large amount of experimental data for Z boson 
decay, which makes reasonable a comparison of bounds set by this novel production process 
with bounds set by other processes. 

The tree-level vertex is absent because it is derived by considering a perturbative 
expansion of the metric about flat space, so that the Z-7-graviton term comes from a per- 
turbation of the Z-7 Lagrangian term, which is zero. (In essence, each vertex involving a 
KK graviton excitation is derived by "hanging" a graviton off a propagator or an existing 



vertex [21, |2^.) Because we are working with bare parameters, we are in effect working 
with a basis for the electroweak sector in which there is Z-7 mixing at one-loop level, but 
this will only enter the Lagrangian via counterterm corrections. It is clear that the ampli- 
tude calculated will not contain terms associated with Z-7-graviton tree-level mixing, since 
we could calculate the amplitude using renormalized parameters instead of counterterms, 
whereby we should have explicitly no Z-7-graviton mixing. 

At leading order, the decay is either into a spin-2 KK excitation or into a spin-0 KK 
excitation (the spin-1 KK excitation does not couple directly to matter [^). The spin- 
2 case is almost identical to the case of Z boson decay to photon plus graviton without 
extra dimensions, which has been considered by Nieves and Pal |l[[. We repeat in this 
paper some of the detail, for sake of completeness. The methodology of the spin-0 case is 
strongly motivated by that of the spin-2 case. 

In many tree-level decays involving particles with small masses, the contribution to 
the decay width of channels involving a spin-0 KK graviton can be ignored, as the ver- 
tices involving the spin-0 KK graviton contain the masses of the other particles (and in 
some cases also contain momentum terms that are zero on-shell). However, the possibility 
of massive particles in the loop means that there is a non-negligible contribution to the 
amplitude in the one-loop calculation from the spin-0 decay channel. 

For the decays into a spin-2 KK graviton and for the decays into a spin-0 KK graviton, 
three types of diagram must be considered. The first two types of diagram are those with 
a fermion in the loop and those with a W boson in the loop. (In principle, the second type 
of diagram also includes diagrams with Goldstone bosons and Fadeev-Popov ghosts in the 
loop, but we adopt the unitary gauge throughout.) A third type of diagram is required 
if we are to work using bare parameters, namely the diagrams containing a counterterm. 
It will turn out that the counterterm diagrams evaluate to zero, which is why we work 
with bare parameters (and not renormalized parameters). The result that these diagrams 
evaluate to zero supports the conclusion in [ffl] that such diagrams need not be considered 
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in the near-flat Standard Model scenario, and supports what appears to be a "miraculous 
cancellation" of divergent terms in the spin-0 KK graviton calculation in this paper. 

We also make use of the argument of Q that, in the spin-2 KK graviton case, con- 
servation of the electromagnetic current and of the energy-momentum tensor implies a 
particular general form taken by the amplitude, which simplifies the calculation. We de- 
rive an analogous argument for the spin-0 KK graviton case. 

It is possible to estimate the form that the decay width will take prior to a full calcu- 
lation. For the decay of a Z boson into a photon and a graviton in 3 + 1-dimensional space, 
Nieves and Pal give the estimate 

r ~ a'^GM^ (1.1) 

by dimensional analysis. They derive this estimate by noting that the graviton coupling 
introduces into the amplitude a factor of k = VSttG, where G is Newton's constant in four 
dimensions; they also note that both the Z coupling and the photon coupling introduce into 
the amplitude a factor of ^/a, where a is the fine-structure constant. (The estimate is then 
obtained by noting that the Z mass, Mz, is the only dimensionful parameter remaining in 
the calculation.) 

For the calculation in the ADD scenario, the factor a^G remains in the width (and 
G is still Newton's constant in four dimensions). However, we sum over the Kaluza-Klein 
excitations of the graviton by using an integral approximation over a "density" p , and 
this density contains a factor i?", where R/2tt is the radius of the extra dimensions and 
n is the number of extra dimensions. The Z mass remains the only other dimensionful 
parameter (since we have summed over all KK graviton masses). We therefore obtain the 
estimate 

r ~ a2(^ii:"M|+™ . (1.2) 
The detailed phenomenology of our more precisely calculated results is covered in 



reference ||12 |. 



The model of toroidal compactification of the extra dimensions with a common com- 
pactification radius R/2tt is something of a toy model, not least because it does not suggest 
a mechanism for confining the Standard Model fields to a 4-brane, as is required for the 
ADD scenario. However, the model still deserves phenomenological investigation, firstly be- 
cause the calculations involved may be useful for understanding models where the topology 
of the compactified dimensions is more complicated, and secondly because the toy model 
still gives some indication of the likely consistency of the ADD scenario with experimental 
observations. 

This paper is organised as follows. In Sections |2|, ^ and ^ we state the Feynman rules 
required and give the diagrams corresponding to a fermion loop, a W boson loop and a 
counterterm, respectively. In Section ^ we present arguments giving the general form of 
the amplitude in each of the spin-2 KK graviton production and the spin-0 KK graviton 
production cases. Consideration of the spin-2 case, identically to that derived by Nieves 
and Pal Q, shows that the amplitude can be expressed in terms of one CP-even and two 
CP-odd coefficients. Consideration of the spin-0 case shows that the amplitude can be 
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expressed in terms of one CP-even and one CP-odd coefficient. In Section we demon- 
strate that the sets of diagrams previously presented satisfy the relevant electromagnetic 
and gravitational transversality conditions, so that it is possible to use the general form 
arguments. We then proceed in Section |^ to calculate the coefficients required for an ex- 
pression of the amplitudes. In Section |8|, we use the expressions for the amplitudes to 
calculate decay widths for the individual KK excitation modes, and then sum over these 
modes (approximating the sum by an integral |^]) to obtain an overall decay width. 

1.1 General notation 

We take Greek indices to range over the four dimensions corresponding to the Standard 
Model brane, and Roman indices to range over the n extra (bulk) dimensions. 
We work with a metric tensor linearised so that 

9\p = r]xp + 2K{hxp + r]xp(l)) (1.3) 

(note that this is a hybrid of the notation of |[l|, |22|, which do not have the extra dimen- 
sional dilaton term, and of [Q], which differs by a factor of two). With this definition, the 
gravitational coupling k, satisfies 

K = V8^, (1.4) 

where G is Newton's constant. The fields generated by the KK reduction may be expanded 
in Fourier modes and redefined in terms of massive fields h^^, A^^ and cp^^ (for the spin- 
2, spin-1 and spin-0 cases respectively) 0. (The spin-1 field does not couple directly to 
matter and is therefore neglected from now on as a higher order contribution.) n is a vector 
giving the excitation level in each of the extra dimensions. 

We note that at each mass level (i.e. for each distinct value of the excitation vector n), 
there are one spin-2 KK graviton excitation h^^ and n — 1 spin-0 KK graviton excitations 
(/)^- to be considered. The n — 1 factor comes from noting that the Standard Model particles 
couple only to the trace (p"^ of the spin-0 particles, and that one degree of freedom is lost 
owing to the linear dependence of the modes 0^- . Equivalently, we may note that vertices 
involving the spin-0 particles always contain a 6ij term, and each external spin-0 particle 
is accompanied by an extra-dimensional "polarisation tensor" Cij, which satisfies the spin 
sum identity 0] 

n{n-l)/2 1 ^ . 1 ^ ^ 



with 



which satisfies 



2 2 

s=l 



riiU 



P^j = - ' (1-6) 



so that when we calculate the modulus-squared of the amplitude, the terms carrying extra 
dimensional indices look like 

(]^Pu'Pfr + \pZ'Pf^ ' (1-8) 
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and, using equation (1/7), this evaluates to n — 1, as expected. 

We define the polarisation tensors £^''{q), e"{k) and £'^{p) as corresponding to the 
spin-2 KK graviton excitation, the photon and the Z boson, respectively. The tensors 
satisfy 

e''{k)ku = Q, (1.9) 

e^(p)p^ = 0, (1.10) 

and 

£'^^\q)qx = Q, £''P{q)qp = 0. (1.11) 
The gravitational polarisation tensor is symmetric and traceless: 

£^%q) = Sf'^q) , £^P{q)r^xp = Q. (1-12) 

The momenta satisfy the on-shell conditions 

p2 = M| , (1.13) 
A;2 = 0, (1.14) 
q^ = ml, (1.15) 

where is the mass of the KK graviton excited to level n, given for the toroidally 
compactified ADD scenario by 

ml = . (1.16) 

Four-momentum conservation {p = h + q) yields the on-shell identity 

2k ■q = Ml -ml. (1.17) 

We introduce notation for the off-shell amplitudes F^J^^{q, k) (in the case of decay to 
a spin-2 KK excitation) and F^{q,k) (in the case of decay to a spin-0 KK excitation), 
defined by 

A^W(g,fe) =f^^*((Z)5-(fe)e^(p)F£) (g,fe) (1.18) 



and 



M^'^\q,k)=e^*{k)e'^{p)Fl£>{q,k), (1.19) 



respectively. 
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2. Fermion loop diagrams 

Figure || contains the diagrams for the Z — > 7 + /i"" process with a fermion in the loop, 
and figure ^ contains the diagrams for the Z ^ 7 + 0" process with a fermion in the loop. 

The Feynman rules for the vertices occurring in the diagrams are given for reference in 
figure ^. Their derivation is given in ^ • Some of the tensor contributions to the vertex 
factors are abbreviated for legibility; the abbreviations used are 

C\pua = r]Xur]pa + VXcrripu - TlXplJua , (2.1) 

cxpf3uih,k2) = {ki ■ k2)Cxppy + r]xpkiyk2i3- 

- [ri\ukipk2i3 + r]Xf3kiuk2p - Vl3ukixk2p + (A ^ p)] , (2.2) 
cfpf^u (ki , k2) = cxppu - MlCxpi3u , (2.3) 

Vxp{kuk2) = \ [-ix{ki + k2)p + 7p(^i + k2)x - 2i]xp{h + h- 2"^/)] • (2-4) 
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We use Qf to denote the fermion charge, and have written 

7^EE7^(Xy + yy75); (2.5) 
Xf and Yf are scalars. We have also defined (as in [^) 



Writing F^pj^l and Ffj,,y to denote the contributions from the fermion loop diagrams to 
the off-shell amplitudes i^^p^^, and respectively, and absorbing constants by defining 

{hf) _ KeQfg (hf) , „s 

^P^"' 2 cos Ow ^^^"^ ' 
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2 cos Owr 




A, 








Figure 3: Feynman rules for vertices in fermion loop diagrams. The extra-dimensional 
indices are omitted from the vertices involving (j)^^ (each vertex gains a coefficient of 5ij). 
Symbols and tensors are defined in the main text. Arrows on bosonic lines indicate direc- 
tions of flow of momenta. 



and 



"' ~ 2cos^vF '''' 



(2.8) 
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we may write the contributions to the amphtudes from individual diagrams as 



Ahf-.a.) 
Xpfiv 

p{hf:h) 
' Xpfiv 

Ahf:c) 

Xpfiv 
p(hf:d) 



(2^)4 



(2^)4 
iaxppa^"u{^) 



Tr [^^S{1 - khuSmxpil, I + q)Sil + q)] , 
Tr [%S{1 - q)Vxp{l, I - q)S{lh,S{l + k)] , 



' Xppu 
p{hf:i) 
' Xpfiu 



Df{k)ctaMmpu{k) 



and 



dH 
dH 



TV 
Tr 



rp{<f>f:c) 

py 
pu 



(2vr)4 
3iu;Upi,{k) , 

rpi'pf.e) _ n 
pv ' 

T(p) = 2a;r?^„M|Df(fc)n;3.(fe), 
where we have written 

1 1 

O-Xppu = ^Xpllpv ~ '^VXp'^pv ~ '^VXu'^pp 



%S{1 - kh,S{l) 0/ + - 2m;) S{1 + q) 
%S{1 - q) - - 2m;) S{lh,S{l + A;) 



and 



Ilpuik) ^ J 



dH 



(2vr)4 

and have defined the propagators 
i(/ + m/) 



Tr[%S{lh,Sil + k)], 



5(1) 



Df{k) 



P 



m. 



-17] 



a/3 



/t2 - M| 



(2.9) 

(2.10) 

(2.11) 
(2.12) 
(2.13) 
(2.14) 

(2.15) 

(2.16) 

(2.17) 
(2.18) 
(2.19) 
(2.20) 

(2.21) 
(2.22) 



(2.23) 

(2.24) 
(2.25) 



for the fermion, photon and Z boson respectively (note that these are different from the 
definitions used in Q). This amounts to a choice of unitary gauge for the electroweak 
sector. (We shall later need the propagator for the W boson, which we take to be that for 
the Z boson with Mz Mw-) 
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3. W loop diagrams 



Figure ^ contains the diagrams for the Z ^ j + h"' process with a W boson in the loop, and 
figure P contains the diagrams for the Z — > 7 + (^" process with a W boson in the loop. The 
Feynman rules for the vertices occurring in these diagrams additional to the fermionic loop 
case are given for reference in figure ^. Their derivation is given in reference 0]. Again, we 
have abbreviated for legibility some of the tensor contributions to the vertex factors; the 
additional abbreviations used are 

Na/S-rih, h, h) = "niS-yih - ^2)0 + 11^a{h - k^)l3 + Tlaisih " ^1)7 , (3.1) 
Cv/3i/(^l) ^2) = Cxpf3u - M^rCxpPu , (3.2) 
d\pal3v{kl, k2, ks) = Cxpa^ikl - k2)u + Cxpauiks - kl)/3 + Cxpf3u{k2 - k3)a + 
+ ['nXa'npu{k2 - k-i)p + r]Xf3llau{k3 - kl)p+ 

+ VXvilapiki - k2) p + {\ ^ p)] , (3.3) 

Ral3pv = '^Tla/sripu " TlapVlSu " 11auVl3p , (3.4) 
Rxpa/Bpu = flXpR-alJup — VXaRp/Bup — VX^Rpaup — VXyRppalS — VXpRpuajJ ■ (3.5) 

Writing F^^J and Fpt^^ to denote the contributions from the W boson loop diagrams to 

the off-shell amplitudes F^^^^ and F^^ , respectively, and this time absorbing constants by 
defining 

f1-'.«. CO, (3.6) 

and 

Fj^'^)^Kegco^ewTli'^\ (3.7) 
we may write the contributions to the amplitudes from individual diagrams as 

X D'§{1 + q)N^s.{l +P,-l- q, -k)D^^{l + p) , (3.8) 
?^1Z^'^ = - / -^,N^Pp{l-k,-l-q,p)D^{l-k)x 

X N^r,{l, -I + k, -k)D'^{l)cY^,{l + q, l)D^{l + q) , (3.9) 

TiZ''^ = i I ■^^N^f.pil,-! -p,p)D^{l)dxparAl +P,-l,-k)D^^{l +p) , (3.10) 

tZ^'^'"^ = i I -^dxpc.pp{l, -I - k,p)Dl^il)N^ru{l + k, -I, -k)D^{l + k) , (3.11) 
TiZ-'^ = -D^\p)cxpsAP,k)x 

X j ^:^N^pp{l,-l-p,p)D-^{l)N,r^{l+p-l-p)D^^{l+p), (3.12) 
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(g) (h) 




(i) (j) 



Figure 4: 1-loop diagrams for the process Z ^ j + K^ involving W bosons in the loop. We 
note that once one has defined a convention for charge and momentum flow in the loop, 
diagrams (a) and (b) must be treated as distinct. 
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Figure 6: Extra Feynman rules for vertices in diagrams with W boson loops. Additional 
vertex functions are defined in the main text. 



X j -^N^sp{l,-l-k,k)D^^{l)N,r,{l + k,-l,-k)Dl^{l + k^ (3.13) 
TiZ'"'' = -^Rap,u j ^Dl{^il)cl,,il + q, + q) , (3.14) 
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^vZ'^ = -icx,au{p. k)D"^{p)R.r,p I ^D^wil) , (3-16) 
TiZ'^ = -^c%,M k)Df{k)R,,p, j j0^D^w{l) , (3.17) 



and 



X N^s,{l+p,-l-q,-k)D^^{l+p), (3.18) 



rj^m--^) ^ _2u;M^^,s I (04^"/3m(^ -k,-l- q,p)D^{l -k)x 



X N^ruil, -I + k, -k)D^^{l)D^{l + q) , (3.19) 



w 

X 



X / (|^^7<5/3(^ - ^> fc)Dj^(/)Ar^,,(/ + k, -I, -k)D^^{l + k) , (3.20) 

Tm-.s) = -2icvM^r,^rRo.p,u I ■j0^D^^{l)D^il + q) , (3.21) 

r(r^^^ = -2iu;Mh,^R^^f,,Df{k) J -^D^^{1) ; (3.22) 

we have omitted a number of lines corresponding to diagrams involving production that 
evaluate to zero. 

4. Counterterm diagrams 

Because we are working at one loop using bare parameters, we must consider corrections 
that arise from a Standard Model Z-7 mixing counterterm. It will turn out that such terms 
give a zero contribution (and this is why we work with bare parameters), and it is sufficient 
just to consider the general form that such terms take in order to demonstrate this. 

Figure |^ contains the diagrams for the Z — > 7 + /i" process that have a counterterm, 
and figure |8| contains the diagrams for the Z — > j+cp"" process that have a counterterm. We 
may divide the diagrams into two classes: those containing a two-point counterterm vertex 
and those containing a three-point counterterm vertex. The two-point counterterm vertex 
is the Z-7 mixing counterterm that occurs in the Standard Model if one uses the tree-level 
diagonalization of the electroweak mixing matrix for one- loop calculations. The three-point 
counterterm vertex arises by considering the gravitational perturbation expansion about 
the Lagrangian term corresponding to the two-point counterterm vertex. 

To derive the Feynman rules for the counterterm vertices in this regime, we need to 
consider the Lagrangian Z-7 counterterm that arises from one-loop renormalization in the 



Standard Model |23|. The relevant terms in the bare Standard Model Lagrangian may be 
written as 

Coz,y = -]^Z^{-r^^^d^ + d^^d')Z, - ^A^i-rj^^'-d^ + d^d^A, + ^^M^Z^t^^^^ Z, , (4.1) 
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(a) (b) 




(c) 



Figure 7: Diagrams for the process Z ^ 7 + /i" involving counterterm vertices. 




(c) 



Figure 8: Diagrams for the process Z — 7 + 0" involving counterterm vertices. 



and applying the renormalization 

zj^ + zj^ Afj, , (4.2) 

A^^Z'J^Z^ + Z'J^A^, (4.3) 

M| ^ M| + (5M| , (4.4) 
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we obtain a mixing counterterm in the Lagrangian, which may be written 



+ (M| + 5Ml) z'ilz'il^^^Z^A- . 



(4.5) 



We may read off from this that the Feynman rule for the two-point Z-7 counterterm vertex 
with momentum k passing through is 



+ 



knki/ 



kfj, ku 



Ml + 8Ml) {zf^Z-Jl) - k- (zyizf, + Z]IIZ]11 



k^ 



[Ml + ml)[z'ilz'il 



We write this contribution in the form 



+ 



(4.6) 



Mk)'ntiu + Bk^ku , 



(4.7) 



where 



Aik) 



[Ml + 6Ml) zT.zTa - (zTzZTa + ZT,ZT^ 



(4.8) 



and 



(4.9) 



For the three-point vertices, we must consider the gravitational coupling expansion of the 
metric. We may write the 0{k) term of the expanded Lagrangian as Q 



(4.10) 



where 



Xp 



9=V 



(4.11) 



and we have replaced Minkowski metric terms 77^1/ in the Lagrangian with the perturbed 
metric Qf^u- We find that the Z-7 mixing terms in the energy- momentum tensor are 



t(Z7) 



'ZZ^ZA 



a/2 7I/2 

'AZ^AA 



X [v^.^mpd'^Z^daA" - r^xpd.Z^^d^A-' - rt^xVupd'^Z^daA''- 
- Vt^pVuxd'^Z'^d^A'^ - VpudxZ^dpA^ - ,]p,dpZ^'dxA^+ 
+r]pxd,Z''dpA'' + r]ppd,Z''dxA'' + r]p,dxZ''dpA'' + r]xudpZ''d^A''] + 

+ (M| + 5M|) Zy^z'J^ [ij^xVpu + VppVXu - Vp.mp] Z^A^ , (4.12) 
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k k 



A{k)r]ai3 + Bk^kp 2ujk {mI + 6M^) Z^J^^^ZA^^^ 
hi 





-K [zj^zjj^ + ZjI^z/j^ cxp^,y{p, k) - Hi (M| + (5M|) zJ^zJ^Cfj_uXp 

Figure 9: Extra Feynman rules for counterterm vertices. Abbreviations are defined in the 
main text. 



with trace 



T\ = -2 (M| + JM|) z'J^z'J^r^p^Z'^A'' 



This yields a Feynman rule for the vertex Z^^{p)A^{k)}i^P''^{q) of 

-K (z'ilz'il + Zilz^lf) cxp,u{p, k) - K (M| + 5Ml) Z'ilz'ilc^ 



(4.13) 



(4.14) 



where C^j^^Xp and c\pp,y{jp,k) are defined in equations ( |2.1| ) and (| 
Feynman rule for the vertex Z'^(p)A'^{k)(p"'{q) of 

2u;k {Ml + 5Ml) Z^^Z^Jlr^p, . 



respectively, and a 



(4.15) 



We note for future reference that the form of the Z-7-/1" vertex is such that there is no 
term in which there are four momenta all carrying Lorentz indices (this will be important 
in showing that the counterterms give no contribution to the amplitude). 
The additional Feynman rules are summarized in figure U. 



With these Feynman rules, and writing F^^^ and F^jt to denote the contributions 

from the counterterm diagrams to the off-shell amplitudes F^^^^ and F^j^ , respectively, we 
may write the contributions from individual diagrams to the amplitudes as 



= '^^p.M k)Df{k) [A{k)r^p, + Bkpk, 



p(hx:h) _ 



z'JzZ'Jl + zTzZ'Il) cxp,u{p, k) + (M| + <5M|) z'J^z'JJC.uXp 

a/3. 



1/2^1/2, 



^ApV. ' = ^CAp/3.(p, k)D''''ip) [A{p)7jp^ + Bp^p 



(4.16) 
(4.17) 
(4.18) 
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and 



F^f = 2ujK7]^^MlDf{k) [A{k)i^fs, + BkfsK] , (4.19) 
F(f ^ _2,^^ (^2 + ^^2) z'J^z'JJv,. , (4.20) 



F(^^^=) = 0. (4.21) 



5. General forms of the amplitudes 

We now derive general forms that must be taken by the amphtudes we are calculating; 
these forms allow simplification of the calculation. The argument for the decay involving a 
spin-2 graviton excitation is that of Nieves and Pal [P; the argument for the decay involving 
a spin-0 graviton excitation is essentially the first part of the argument for the spin-2 case, 
and is given first as it is more straightforward. 

The arguments rely upon Ward-Takahashi identities that are consequences of electro- 
magnetic and gravitational gauge invariance. This gauge invariance is shown in the next 
section for each of the sets of fermion loop diagrams, W loop diagrams and counterterm 
diagrams separately. 

5.1 Decay into spin-0 excitation and photon 

Electromagnetic gauge invariance implies conservation of the electromagnetic current, which 
(transforming to momentum space) gives 

k''F(t\q,k)=0. (5.1) 
We may expand F^^^ about A; = 0, writing 

with T^^ independent of k. Equation ( |5.1| ) then implies that 

fe%'l = 0, (5.3) 
^"^%U = 0, (5.4) 

and since this is true for all k with I/cqI ^ in the centre of mass frame and /c^ = 0, we 
can deduce that 



and 



Tl = (5.5) 



This means we may write the on-shell amplitude ( 1.19| ) as 



M^'f'\q,k)=e^Jp)n%,o., (5.7) 
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where 



and 



Recalling equations ( |1.9[ ) and ( 1.14| ) {e^{k)ky = and A;^ = 0), along with the momentum 



conservation relation p = k + q, and considering terms that could be contained in t^i/Q., we 
see that the only terms contributing to the amplitude are ones not involving a Levi-Civita 
tensor, of the form 

scalar x 4(p)(A;'^e"* - A;"e''*)(r?^^(?„ - rj^^q^) , (5.10) 

and ones involving a Levi-Civita tensor, of the form 

scalar x ef^{p)e''°''^\k^es* - ks£^^){r]^uqa - -q^aqv) ■ (5.11) 

Relabelling indices so that the polarisation tensors may be written as common coefficients 
of the overall amplitude, we may rewrite the above as 



Flf = (k^Qu - k ■ gr?^.)FW + (e^.^^g"^^)^^'^) , (5.12 



where F^'^^ and are Lorentz scalar s. 

5.2 Decay into spin-2 excitation and photon 

The argument in the case of decay involving a spin-2 graviton excitation may be viewed 
as an extension of the case involving a spin-0 graviton excitation. It is given in full in 
reference IQ for the case of a massless graviton, and is sketched here. 
Similarly to the spin-0 case, we may write 

Fxpi_iu = "^Xp^iv + k Ty^pp^a 1 (5.13) 

and we may use the condition k^F\ppy to deduce that we may write the on-shell amplitude 



( 1.18 ) in the form 

A^W(g, A;) = e^P{q)e^^{p)n*txpp,ua , (5.14) 



where f^a is as defined in equation (|5.8D . 

We may also use gravitational gauge invariance to derive a manner of writing the am- 
plitude, in a similar fashion to the way we have used electromagnetic invariance (although 
the details are complicated by the extra Lorentz index in the gravitational case). 

Writing 

M^''^=£^''{q)jxp{q,k), (5.15) 

and expanding 

jXp = fxp + Q^jlpa + q^q^jlpar , (5.16) 
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where j^^ and j^^^ are independent of we may use the conditions 

q^jxp = 0, q'jxp = (5.17) 
(following from conservation of the energy-momentum tensor) to deduce 

Jap = 0, (5.18) 
jipa = 0. (5.19) 

In addition, has a number of symmetries arising from its definition (namely \ ^ p 

symmetry, a ^ t symmetry and antisymmetry under interchange of either of (A, p) with 
either of (a, r)). We also note that no term in the expansion of j\p„^ contains a qx or qp 
term, because such terms vanish on contraction with the gravitational polarisation tensor 
(equation (|1.11| )). Further, we need not consider terms in q^ or g,-, since q^ = mi, and 
so such terms reduce to lower order terms in equation ( 5.16| ). Combining these properties 



with the expression for the form of the amplitude obtained from consideration of the 
electromagnetic gauge invariance yields a general form of the amplitude that may be written 
as 

^ipi. ={ikxqu - k ■ qvux){kpq^, - k ■ qv^.p)F^''^ + exuaf^q^'k^ikpq^, - k ■ W^p)Ff ) + 

+ {kxqu - k ■ gr?.A)epM"/39°^^^i''^ } + (A ^ p) . (5.20) 



6. Transversality conditions 

In order to apply to any set of diagrams the arguments presented in section |5| regarding the 
general forms of the amplitudes, we must establish that if we take the contribution given 
by the set of diagrams to the amplitude and contract it with the photon momentum, this 
yields zero. (In the case of the decays involving a spin-2 graviton excitation, we also need 
to establish that contraction with the graviton momentum yields zero.) 

6.1 Electromagnetic transversality: spin-2 case 

We begin by considering the case of decays involving production of a spin-2 graviton ex- 
citation. In the spin-2 case, the contributions to the amplitude of the diagrams involving 
a fermion loop and of the diagrams involving a W boson loop are essentially the same as 
the contributions given in [||]. We repeat here for completeness an outline of the derivation 
of the transversality results in those cases. The counterterm diagrams are not considered 
explicitly in Q, and we present more fully the argument for these diagrams. 

6.1.1 Fermion loop diagrams 

Using the identity^ 

k''S{l)j^S{l + k)= iS{l) - iS{l + k) (6.1) 

^The identities differ from those given in M by factors of i owing to the different propagator definitions. 
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with equation ( |2.22| ) implies that 

rn^^(A:)=0, (6.2) 
from which, using equations ( 2.12| ) and ( |2.14D , it fohows that 



^^r£i? = 0, k^T^ = ^. (6.3) 



Using the definition of equation (|2.2| ) , we also have 

A:"cAp/3.(p,/c) = 0, (6.4) 

from which it follows that 

rTlJ^-) = 0. (6.5) 

We can show that 

^^W.^ + k'^T^^/^^ + k^T^^ = , (6.6) 

by using equation ( |6.1D on the sum k'^T^xJ^j^v' + ^'^'^Xpiu^ ^^'^ writing the resultant term so 
that the only propagator terms are S{1) and S{1 + p), and then using the Ward-Takahashi 
identity 

Vxp{l + q,l)- Vxpil + p,l + k) = axpapl^'k" . (6.7) 



This establishes that 



k^Ti%l = 0. (6.8) 



6.1.2 W loop diagrams 

Electromagnetic transversality of the W loop diagrams contributing to the spin-2 excitation 
producing decay may be demonstrated by using a number of Ward-Takahashi identities. 
The identities 



r^D'^il + r)N^(3-yil + r, -I, -r)D^{l) = iD^{l) - iD^il + r) , (6.1 
cYpafiil + - c^pap{l + r + s, / + r) = r'^dxpafSyil + r + s, -I, -r) , 



.10) 



and 



Na,f3,'yi^ - r,-l - S,r + s) - Na,p,'y{l, -I - r - S,r + s) = Ra/S-yS , (6-11) 

together with equation (6^), can be used with equations ( ^^ ) to ( |3.17| ) to deduce that 



l,urp{hW:a) ,urp{hW:h) rurp{hW:c) i.urp{hW:g) _ p, 1 9^ 

^Xppu +'^^Xppu ^'^^Xppu '^'^^Xppu [G.ll) 

k^TiZ-'^ + k^Til^-^^ = 0, (6.13) 

k'TiZ'^ = 0, (6.14) 

k^TiZ'^ + k^Ti:;^-^^ = 0, (6.15) 
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and 



and therefore that 



6.1.3 Counterterm diagrams 

Using equation (|6.4D with the definition of equation ([4.18|) impHes that 



k''Tt^J'^ = Q, (6.16) 



^^^^£^ = 0. (6.17) 



™o uo o.x.p..^ ... ....... u. . .^('^""^ -'^ 



Equation ( |6.4[ ) also allows us to simplify the forms of Fj^p^y and Fj^p^y ■ If we expand 



the definitions of cf (p,k) and D^^ {k) from equations ( ^ ) and ( p. 25 ) respectively, this 



yields 

= '^^(^)(^^^Ap - k^r^,p - kp^p,) = -FF^J;;^) , (6.19) 
from which we deduce that 

6.2 Gravitational transversality: spin-2 case 
6.2.1 Fermion loop diagrams 

Gravitational transversality of the fermion loop diagrams contributing to the spin-2 gravi- 
ton excitation production process may be seen through use of the identities 

S{l)q^Vxp{l + q, l)S{l + q)= i{lp + qp)S{l + q) - ilpS{l) + 

+ \ (Tp^ - ilp) S{1 + q)- '-S{1) (7p^ - ^7p) (6.21) 

and 



7a7/37p + 7p7/37a = 2 (r?a/37p + f?/3p7a " f?ap7/3) , (6-22) 

together with the momentum conditions and with changes in the integration variable in 
some of the loop momenta integrations. Fuller details may be found in reference [|l|. 

6.2.2 W loop diagrams 

If we use cyclic identity 

Rap-yS + RfS-yaS + R'ja/SS = , (6.23) 

together with the identities 

Nap^{l, -I - s,s) - Naf3^{l, -I - s + r,s-r) = r^Rp^aS (6-24) 
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and 



-I - s,r + s) = r^R 



(6.25) 



the identity of equation ( 5.10 ) and the identity 



q^D^il + q)cY^^pD^^ = i{lp + qp)D^{l + q) - UpD^il)- 



zr^^{q'^D^^{l)+q^D'g{l + q)) 



then it is possible to derive the result 



0. 



(6.26) 



(6.27) 



More details of the manipulations involved may be found in reference 
6.2.3 Counterterm diagrams 

Gravitational transversality of the counterterm diagrams may be shown simply by con- 
tracting each of the diagrams with q^ and summing the three resultant expressions, if one 
uses the on-shell momentum conditions for the Z boson and the photon (equations ( 1.13| ) 
and (|1|)) H. 

6.3 Electromagnetic transversality: spin-0 case 

We now turn to demonstrating electromagnetic transversality of the separate classes of 
diagrams in the case of decays involving spin-0 graviton excitation production. This is 
more straightforward than the spin-2 case, because there are fewer non-zero diagrams to 
consider. (In addition, we only need to demonstrate electromagnetic transversality to 
establish that the contributions from the classes of diagrams are of the general form given 
in equation ( |5.12| ): there is no corresponding gravitational transversality condition to be 
satisfied). We use some of the identities from the spin-2 case. 



6.3.1 Fermion loop diagrams 

Using equation ( |6.2D with equations ( p. 18 ) and ( |2.20| ) implies that 



- Ill/ 



0. 



In addition, we can use equation ( |6.1D with equations ( 2.15 ) and ( |2.16| ) to write 

dH 



and 



(2vr)' 



dH 



rTr 



% {3(1 -k)- 5(0} [b + -J- 2m f ) S{1 + q) 



(6.28) 



(6.29) 



Tr 



%S{1 - g) ( 2/ 



2m/ {S{1)- S{l + k)] 



(6.30) 
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and we may redefine the loop momenta in the integrands to obtain 

dH 



so that overall 



(6.31) 
(6.32) 

(6.33) 



6.3.2 W loop diagrams 

We begin by taking the expression for fe-T^r^'^ using equation ( |3.20D , and applying 



dH 



equation (|6.9| ), obtaining 



(2vr)' 



N^Sfiil, -I - k, k) \Dt{l + k)- D^^il)] . (6.34) 



We then redefine the loop momentum integration variable so that the integrand has a 
common coefficient L'^(/), and apply equation (|6.11| ), to obtain 



kurj^m-.f) ^ 2iu:Mlrj^^k-R^si5uDf{k) 



dH 



(2vr) 



(6.35) 
(6.36) 



We also redefine the loop momenta integration variables in the expressions for 



{4>W:a.) 



and 



(equations ( 3.18|) and (|3.19|) , taking I ^ I — q and I ^ l + k, respectively), and use 



equation (|6.9|) , to write 



D^^{l + k)-D'^{l) 



Pa, 



+ 



D<^{l + k)-D^{l) 



2iujMu 



dH 

(2vr) 



^D^{1) [N^(3f,{l -k,-l + k-p,py 



-N^f3f,{l,-l-p,p)]D^^{l + q) (6.38) 



(6.37) 



dH 
(2^ 



D^{l)D^{l + q) 



= 2iujM^^r]^rk''Ral3fiiy 

~ ^lU 1 

where to obtain equation (|6.39| ) we have used equation (|6.11| ) . 

As the other diagrams evaluate to zero, we are now able to establish the result 



(6.39) 
(6.40) 

(6.41) 
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6.3.3 Counterterm diagrams 

Contracting with k'^ the expressions for pjit^'^^ and Pj^^'^^ given in equations ( |4.1S| ) 
and ( [4.20I ) respectively, and recaUing the expressions for A{k) and B given in equations ( |4. 
and (H^) respectively, along with the expression for the propag ator Df{k) given m equa- 
tion (|2.25|) , it is straightforward to see that 

Given that F^^""^'^ = 0, we can therefore deduce that 

rF^^x) = 0. (6.43) 

7. Calculation of the amplitude coefficients 

Having shown that the individual sets of diagrams (ones with fermion loops, ones with 
W boson loops, and ones with counterterms) separately satisfy the requisite conditions 
to give a contribution to the amplitudes matching the given general forms, we can now 
calculate these contributions separately. We begin by showing that the contribution to 
the amplitude from the counterterm diagrams is zero in both the decay involving a spin-2 
graviton excitation and the decay involving a spin-0 graviton excitation. This will mean 
that the contributions from the other sets of diagrams must be finite, and so we should not 
be surprised when we see a "miraculous" cancellation of divergences in the case of decay 
to a spin-0 graviton excitation. Indeed, we note that as we have not made any particular 
assumptions about the number of fermions, we know that there cannot be a cancellation of 
infinities between the fermion loop diagrams and the W boson loop diagrams. This means 
that the contribution from each of the individual sets of diagrams must be finite. 

7.1 Counterterm diagrams contribution 
7.1.1 Spin-2 case 

We note that there are no Levi-Civita tensors in any of the counterterm-containing dia- 
grams that contribute to the decay amplitude in the spin-2 graviton excitation case (equa- 



tions (|1|) to (|4l8|) ). The dia grams therefore give no contribution to either or Fr 



2 



and we need only consider the contribution to F^'^\ 

To consider the contribution to F^'^^ we need consider only one of the terms of which 
it is a coefficient, and (as in ||l|) we choose the term F^^^kxkpq^^qp (noting that this term 
appears twice in the expression given in equation ( |5.20| )). Equation ( |1.10| ) implies that 

^zW^M = -ezb)^M> (7-1) 

so that terms requiring consideration also arise from terms of the form k\kpkfj_qij. 

For diagram (a), we note that the contribution has an overall coefficient of k^, and 
therefore does not contribute to the amplitude. For diagram (b), we note that there is no 
term containing four momenta with Lorentz indices, so there is no contribution from this 
term either. Similarly for diagram (c), once we note that the term with a B coefficient 
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contains a factor of so provides no contribution to the amplitude, we can see that there 
is no remaining term with four Lorentz index-carrying momenta. The contribution to F^'^^ 
from the counterterm diagrams is therefore zero. 

7.1.2 Spin-0 case 

We note that, as in the spin-2 case, there are no Levi-Civita tensors in any of the counterterm- 
containing diagrams that contribute to the decay amphtude in the spin-0 graviton excita- 
tion case (equations ( [4.19 ) to ( 4.21D ). The diagrams therefore give no contribution to 



and we need only consider the contribution to F^'^^ 

Similarly to the spin-2 case, we need consider only one of the terms of which F^'^'^ is 
a coefficient, and we choose the term k^j^q^. In this case, the term is not repeated in the 
expression for the general form of the amplitude. Again, equation ( 1.10| ) implies that terms 
requiring consideration also arise from terms of the form Q/j^qu- 

It is sufficient for the counterterm diagrams to note that no diagram expression depends 
upon the momentum q of the spin-0 particle, so that there is no contribution from the 
counterterms to the amplitude. 

7.2 Loop diagram contributions to the spin-2 amplitude coefficients 

We shall see that there are no Levi-Civita terms to be considered, so that the contributions 
to Fi^^ and ) are zero. In considering contributions to F^^^ , we shall again look at terms 
of the form kxkpq^q^ (or kxkpk^q^). 

7.2.1 Fermion loop diagrams 

The determination of the contribution from diagrams containing fermion loops is very 
similar to the massless case considered in 

With respect to the contributions to f'^'''^ and F2^\ we note that for diagrams (c), (d), 
(e) and (f), any term with a Levi-Civita tensor also has a metric tensor symmetric in two 
of the indices of the Levi-Civita tensor, so that there is no contribution from any of these 
diagrams. 

With respect to the contributions to F'^^\ we note that none of diagrams (c), (d), (e) 
and (f ) contributes a term of the relevant form: this is straightforward to see for diagrams 
(c) and (d), which only depend upon p and k respectively; it is relatively easy to see 
for diagrams (e) and (f) if one uses the definitions of c\ppy{p,k) and cfp^^{p,k) given in 
equations ( |2.2| ) and (^), respectively. 

It is therefore necessary only to consider contributions from diagrams (a) and (b). The 
contributions from diagrams (a) and (b) may be written in the form 



rp{hf:a) 



rp{hf:h) 
Xpfiu 



(2vr)^ 



(/ - ky -mj {1 + qY - mj (^P - mj^ 



dH IpTr [%(/ - ^ + m/)7A(/ + mj)7^(/ + ^ + mf)] 



(2vr) 



(/ + ky - mj {I - qY - mj [l"^ - mj^ 



(7.2) 
(7.3) 
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As in reference iQ we can use the cyclic property of the trace and charge conjugation 
relations 

r'-i^,M^,5r:< _ (^m^,^\'^ (7.4) 



to write the sum of these terms as 



where 



(7.5) 



(7.6) 



As the trace does not contain a 75 term, the contribution to the amplitude will not contain 



a Levi-Civita tensor, and so the contributions to F^'' and F2''' are zero. A Feynman 
parameterization yields 



Ah) 



rp{hf:a) rp{hf:h) 
Xppu ' Xppu 



dH 
(2^ 



dx 



dy 



fxppuil + xk-yq) 



m?f + y{l 



X — y)m'i + xyM\ 



(7.7) 



and considering only the terms we have previously noted, there are no divergences to be 
considered and we can integrate with respect to the loop momentum, obtaining a contri- 
bution to F^^^ from each fermion in the theory of 



p{hf) 



Keg 



where 



J{X,Y,Z) 



cos 9w 

1 j-l—X 

dx / dy 
Jo 



QfXfJ{mf,mfi,Mz) . 



x'^y{l 



y) 



X^-y{l 



y)Y^ — xyZ^ 



(7.8) 



(7.9) 



7.2.2 W loop diagrams 

The determination of contributions from diagrams containing W boson loops is also very 
similar to the massless case considered in 

There are no Levi-Civita terms in the individual diagram contributions, so the overall 
contributions to F^'^^ and -Fg'*^ are zero. We need therefore only consider the k\kpqfj_q^-like 
terms to obtain a contribution to the coefficient F^'^\ 

The contributions from diagrams (g) and (h) have no dependence on k, so cannot take 
the form we are considering, and therefore may be ignored. Similarly to the fermion case, 



the forms of c\pav{p,k) and c^p^^(P5^), given in equations ( |2.2D and ( |2.3D , respectively, 
are such that the combinations T^^^^^ +^Ap)h^''* ^"^^ '^Xp^u^^ '^'^Xp^i/^^ contain no relevant 
contributions to the coefficient F^'^^ (see Q for more details). Expanding the expressions 
for diagrams (c) and (d), we see that no terms in the expressions contain a term, so these 
diagrams provide no relevant contribution. We are therefore left with only the contributions 
from diagrams (a) and (b). The contributions from these diagrams are equal: this may 
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be seen by applying a change of variable in one of the integrands in the expressions for 
'^Xp^u^^ and T^^^i^'^\ or perhaps more straightforwardly by observing that the contributions 
obtained from the diagrams are independent of the direction in which the W boson loop is 
traversed. 

It may be shown after some manipulation Q that the contribution from diagram (b) 
is contained in the equation 

-M^-^-y) ^ 

(2vr)4 Jq 7o [/2 - + y{l-x- y)ml + xyM^] ^ 

+ non-contributing terms , (7-10) 

and therefore, after performing the integration over the loop momentum, that the overall 
contribution from the W boson loop diagrams to the coefficient F^'^^ is 

F'""' = ^^^^ [s - ^) J(M„„m.,. M.) . (7.11) 

where J{X,Y, Z) is defined in equation (|7.9|) . 



7.3 Overall contribution to the spin-2 amplitude coefficients 

We have determined so far that the coefficients F^^^ and F2^^ are zero. We have determined 
also that the coefficient F^'^^ is given by 



Keg 



Att"^ cos 9\Y 

cos^Ow ^og2g^ ) J{Mw,m^,Mz) - 2J2Q fXfJ{mf,m^, Mz) 

■(7.12) 

To proceed, we approximate the analytically intractable integrals of the form J{X, Y, Z). 
We note that all the integrals in which we are interested take the form J{X, nifi, Mz), and 
since the mass of the Kaluza-Klein mode is constrained by < Mz, we may consider the 
integral only in the cases where <Y < Z . 

We consider the expression for J{X,Y,Z) given in equation ([7.9D, and repeated here 
for convenience: 

J(X. Y. Z) = /' ,,__0i^pL_ . (7.13) 







We note that over the range of integration, the maximum value of the expression xy is 
1/4, and the maximum value of the expression y(l — x — y) is also 1/4. This allows us 
to consider two limiting cases for approximation of the integral: X/iY + Z) ^ 1/4 and 
X/{Y + Z) <C 1/4. With the constraint < y < Z, it is sufficient to consider the cases 
X/Z > 1/2 and X/Z < 1/4. 
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To evaluate the case XjZ ^ 1/2, we may approximate the integral by J(X, 0, 0), which 
is relatively straightforward to calculate, and yields the result 



J(X,0,0) 



(7.14) 



360X2 • 

To evaluate the case XjZ <^ 1/2, we may approximate the integral by J(0, Y, Z). We 
may change the order of integration so as to perform the x integral first; this is possible 
analytically. We obtain an integral polynomial in (1 — y) that we can evaluate analytically. 
Evaluating this integral, we obtain 

1 Z^ Y'^Z'' Y^Z^ 

j(o,y,z) 



■ + ■ 



12(Z2 - y2) 8(^2 _ y2)2 4(^2 _ y2)3 4(^2 _ y2)4 
To proceed, we use the identity 

^2 ^ C2 _ 

C^-B^ = C2 - B2 ~ 
to eliminate all terms with a numerator, as well as the identity 

C2-52^ 



log 



/y2 
1^2 



log 1 



C2 



(7.15) 



(7.16) 



(7.17) 



and expand the logarithm term (the expansion is valid for < Y < Z , and the final answer 
turns out to be valid for y = as well). We obtain 



1 °° 1 
^^°'^'^) = -2g(, + l)(,+2)(,+3) 



z-2i(z2-y2)i-i. 



(7.18) 



We now return to our expression for F^^^ given in equation ( |7.12| ). We approximate 
the integrals by using the form of equation (7.14) for J{M\y,mfi, Mz) and J{mi, rufi, Mz), 
and the form of equation ( 7.18| ) for the integrals relating to fermions other than the top 
quark. This gives 



Fih) 



Keg 



47r2 cos 6'vi/M2 



1 

360 



1 



cos2 



+ ( 5 - ^ sin^ 9w 




(j + 2)(j + 3)(j + 4) 



Ml 



Ml 



1 



180 135 



sin^ 9w 



(7.19) 



(where we have used that M\y = Mz cosOw)- Using e = ^sin^iy and sin^ 9w = 0.23, 
along with a numerical value for Mz /mt p5| , we obtain the approximation 



Ke 



47r2M2 



0.030 + 4.6 I 

j=0 



1 



2 \ J 



1 

(j + 2)(j + 3)(i + 4) V 



(7.20) 



This approximation agrees with that given in reference ||T| in the case = 0. 
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7.4 Loop diagram contributions to the spin-0 amplitude coefficients 



As in the spin-2 case, we shall see that there are no Levi-Civita terms to be considered, 
so that the contributions to F^"^^ are zero. In considering contributions to F^'^\ we shall 
again look at terms of the form k^Qi,. 

7.4.1 Fermion loop diagrams 

The contributions from diagrams (d) and (f) are zero, since the contributions only depend 
upon k, so they are not of a form to contribute to F^'^^ and symmetry considerations show 
that they do not contribute to Fi'^\ A similar argument shows that diagram (c), which 
depends only upon p, does not contribute to either of the coefficients. 

As diagram (e) trivially gives a zero contribution, it follows that it is necessary only 
to consider contributions from diagrams (a) and (b). 

The methodology employed is very similar to that of the spin-2 case, developed in Q. 
We begin by writing out expressions for the fermion propagators, so that we may write the 
contributions as 



rp{4>f:a) 



-2iuj 



-2iuj 



dH TV (/ - ^ + mj) 7^ (/ + mf) (§/ + f ^ - 2m/) (/ + ^ + mj)] 



{l + qf 



m. 



{I - kf 



m. 



m. 



(7.21) 

dH TV (/ - ^ + m/) (|/ -y- 2mf) {t + mf)^^{t + $ + m/)] 



{2^Y 



{l + kf 



(/ - qf 



rrij 



(7.22) 



For the second integral, we change the integration parameter to — /, and apply the cyclic 
property of the trace, the trace-reversal invariance of strings of gamma matrices, the in- 
variance of transposing matrices inside the trace, in addition to the charge conjugation 
properties of the matrices. We are able to obtain a term in Yj that cancels with the cor- 
responding term in the first integral (so that there is no Levi-Civita term and therefore 
no contribution to F^"^^), and a term in Xj equal to the corresponding term in the first 
integral. We are therefore able to write 



dH 



{2nY 



fUl) 



{1 - kf 



{l + qf 



m. 



(7.23) 



where 



V(0 = IV 



7m (/ + f 



rrii 



3 , 3 
2^+4^ 



2m/ ) (/ + mf) (/ 



f. -I- m/) 

(7.24) 
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Applying a Feynman parameterization and continuing to D dimensions, the equation be- 
comes 



flU ' 111/ 



(2vr) 



D 



JO 



P — m'j + y(l — X — y)m'i + xyM"^ 



(7.25) 



We are interested in terms that wih contribute to the coefficient of A;^g^. We may use the 
symmetry properties of the / integral to discard terms that are odd in /, and to replace 

lah J^Valsl"^ ) (7.26) 



(see for example page 477 of reference [^6[). We find [^] that the term in the numerator 
of the integral containing k^q^ may be written as 



SMzxy (2 - 3x - 6y + Axy + 4y^) + Smly {1 - x - y) {l - 3x - Ay + 4xy + 4y^) + 
+ mj{5-3x- 20y + 20xy + 20y^) + 



+3/2 (^{-1 + x + 4y-4xy- 4y^) + (-1 + 3x + 8y - 8xy - Sy^)^ 



(7.28) 



We note that the integral appears to have a divergent term. However, it will turn out that 
the coefficient of the divergent term is zero after integration over the Feynman parameters. 

It is possible, although not entirely straightforward, to evaluate the integral with the 
numerator given above. In order to do so, we adopt the following strategy: 

• Perform the integrations over the loop momentum /. For the numerator terms not 
containing factors of /, we use the standard result. For the numerator terms con- 
taining factors of /, we use the results given in Appendix^. We obtain a result that 
appears to diverge as D — > 4. 

• Order the Feynman parameter integration so as to perform the x integral first. 

• Perform the x and y integration for the coefficients of the terms that appear to be 
D ^ 4 divergent. (These are also the coefficients of 7 and log(47r) "MS-like" terms, 
which we therefore consider at the same time.) The terms vanish. 

• To simplify the algebra, rewrite the numerator of the l/lm'j, — y{l — x — y)m'i — xyM'^] 
term to eliminate terms with numerator coefficients involving m'j. (Note that doing 
this before approximating the integral in the way done below makes no difference 
when making the light fermion approximation given, and is in fact equivalent to 
including an extra term in the series expansion for the fraction when making the 
heavy fermion approximation given.) 



-32- 



At this stage, there are terms in the integrand that do not contain masses (these 
arise from the step above, and from the "anomalous" terms arising from the loop 
integration that looks like equation (|A.1[ )). Integrate these terms with respect to the 
Feynman parameters. 



Integrate the logarithmic term with respect to x by parts, obtaining a contribution 
to the l/[mj — y(l — x — y)m'i — xyAl'^] term and a logarithmic term to be integrated 
with respect to y only. 



Integrate the remaining logarithmic term with respect to y by parts. One of the re- 
sultant terms is zero. Approximate the denominator of the other term by considering 
the cases <C M'^/A (in which case take mj = 0) and -^1/4 (in which case 

take M| = 0). 



In order to perform the remaining integral (the integrand has a denominator of 



m 



f 



2/(1 



y)rni — xyM^]), approximate the denominator of the integrand by 



considering the cases mj <C (in which case take mj = 0) and m j ^ (in 

which case take = 0, = 0). The latter approximation leads to a relatively 
straightforward integral. For the former approximation, the integral may be evalu- 
ated by considering the numerator and denominator as polynomials in x, and writing 
the integrand in quotient -|-remainder form. 



Following this strategy, we obtain that for each fermion, the contribution to the amplitude 
coefficient F(<^) IS 







8 I ^ 
3 ^ 45 



7m?.+llM| 



for mj < iM| , 
for » iM| . 



(7.29) 



7.4.2 W loop diagrams 

There are no Levi-Civita terms in the individual diagram contributions, so the overall 
contribution to is zero. We need therefore only consider the k^qiy-like terms to obtain 
a contribution to the coefficient F^'l'\ 

The contributions from diagrams (f) and (j) are zero, since the contributions only 
depend upon k, and therefore do not yield a term of the form we are considering. 

With the exception of diagrams (a), (b) and (g), the other diagrams trivially give a 
zero contribution. We need to consider the contributions from each of diagrams (a), (b) 
and (g). 

We begin by considering the contributions from diagrams (a) and (b). These diagrams 
yield an identical contribution. We may see this through consideration of the integral 
that contributes to T^^^''^\ given in equation ( |3.19| ). If we take this integral, change the 
integration parameter taking I ^ —I — q, relabel the dummy suffices and note that 



<(0 = ^^"H) 



(7.30) 
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and 

Naf)j{ki,k2,h) = -Npaj{k2,ki,k3) , (7.31) 

we are able to obtain the expression for tI^^'-^\ We can therefore obtain the contribution 
from both diagrams by considering the form of the contribution from one of the diagrams, 
and we choose diagram (b). We write the contribution in the form 

M^J (27r)4 [(I - kf - M^] [{I + qf - M^] [P - M^] ' ^'-'"^ 

where 

Sf,^(J,) = My^r}^sNai3fji{l -k,-l- q,p)N„ruil, -I + k, -k)x 



(7.33) 

Applying a Feynman parameterization and continuing to D dimensions, the equation be- 
comes 



J (27r)^7o Jo 



X . . , w., S^^il + xk-yq) 



[/2 _ + y(l - X - y)m| + xyMl] ' ' 

(7.34) 

We leave the equation in this form for now, and turn to the contribution from diagram 
(g). Unlike the corresponding diagram for spin-2 particle decay, this diagram provides a 
contribution to the amplitude coefficient. We write this contribution in the form 

M^[P-M^^][il + qr-M^]' ^'-''^ 

where 

We wish to combine this contribution with that for diagrams (a) and (b), so we write it in 
the form 

rm-.,) f u,Ai) [il-k)'-M^] 



J (27r)4 [P _ M^] [{1 + g)2 - M^] [{I - kf - Ml,] ' 

Applying a Feynman parameterization and continuing to D dimensions, the equation be- 
comes 



J WW Jo ''Jo 



i-x U^uil + xk-yq) (l+{x-l)k-yqy -M^ 



[P -M^ + y{l-x- y)ml + xyM|] ^ 

(7.38) 
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We may now combine the contributions from diagrams (a), (b) and (g). Using the symmetry 
properties of the / integral given in equations ( 7.26D and ( |7.27D , and taking a factor of 
Aiu/My^ outside the integral, we find |24] that the term in the numerator of the integrand 
containing kfj^q^ may be written as 



M| ( 2xy2 - 2xy^ - ^x'^y - x'^y^ + ^x^y - x^y' 



+ mlMl [ 2y^ - Ay" + ly" - 3xy^ + Sxy^ + x^y 
-^x^y^ + 2x'^y^ - x^y + 2x^y'^] + 



xy^ - ^x^y + ^x^y^ - 2x^y^ + ^x^y - x^y^ ) + 
1 



+ M^M| ^-2 + 4y- 6y^ -2x + 2xy + -x^ - bx^y ) + 

+ My^rtni ^-5xy + 5xy^ - + bx^y 
+ (12 + 8y - 8y2 _ 4^ - 8xy + (-8y + Sy^ + 8xy)) + 



+ lH Mi 



-2 + 2y^ + 2x + xy - -x^ + 2x^y ) + 



+ — (2 - 4y + 4y^ + X + xy - x^ 



+ 4x2y)^ + 



+ m. 



2xy - 2xy^ + -x^ - 2x2y ) + 



+ 



D 



-y + — X + 5xy — 4xy^ + x^ — 4x^y) ) + 



+ (4 + 5x) + - (-12 + 5x) 



D 



-x+ — (1 - 2x) 



(7.39) 



We evaluate this integral using a similar strategy to that for the fermion loop case. To 
simplify the algebra, we rewrite the numerators of some fractional terms in order to ob- 
tain fewer terms requiring the application of approximations. Subdominant terms will be 
affected by the order in which rewriting a numerator and approximating the denominator 
is performed. However, this effect will be negligible at the order of the approximations 
made to the fractions themselves. Nevertheless, we record the order in which we proceed 
through the algebra, so that it is possible to reproduce the result. 
The strategy adopted is as follows: 



Perform the integrations over the loop momentum /. As in the fermion case, for the 
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numerator terms not containing factors of we use the standard result, and for the 
numerator terms containing factors of Z, we use the results given in Appendix 
Again we obtain a result that, at first, appears to diverge as — > 4. 

• Order the Feynman parameter integration so as to perform the x integral first. 

• Perform the x and y integration for the coefficients of the terms that appear to be 
D ^ 4: divergent (and of the terms that are "MS-like"). The terms vanish. 

• Rewrite the numerator of the 1/ [M^/ — y{l — x — y)m'i — xyM^'j term to eliminate 
terms with numerator coefficients involving M^^ (as in the fermionic case, when taken 
with the approximation for the denominator below this is equivalent to adding higher 
order corrections in the series expansion for the fraction). 

• Integrate the terms that are non-fractional and do not contain a logarithm. 

• Integrate the logarithmic term with respect to x by parts, obtaining a contribution to 
the 1/ — y(l — x — y)?n.^ — xyM^] term and a logarithmic term to be integrated 
with respect to y only. Rewrite the numerator of the fraction to eliminate terms with 
coefficients involving M^r. 

• Integrate the remaining logarithmic term with respect to y by parts. The only term 
that results is an integral with fractional integrand. This term can be integrated if we 
approximate the denominator by M^r■ (Note that the result below does not involve 

w 



rewriting the numerator term with coefficient for this term only. 



• Perform the remaining double integral by approximating the denominator of the 
fractional integrand by M^^. 

Following this strategy, we find that the contribution from the W boson loop diagrams to 
the coefficient F^'l''^ is 



= Keg cos Ow 



UJ 



14 + iiZ^ _ + i^Zi . 1 _ 3 M|l 

15 15 210 210 20 



(7.40) 



We have determined so far that the coefficient F^'^'^ is zero. We have also determined that 



7.5 Overall contribution to the spin-0 amplitude coefficients 

We have determined so far tha 
the coefficient F^''^^ is given by 

1 7ml + llMl 



F^"^) = Keg cos Ow 



(47r)2 
KeQtg ioXt 



11 ml 37 M| 19 mi 1 m|M| 3 M| 



14 J n ^ J !3_ + 

15 15 210 210 20 



2 cos 9w (47r)^ 



3 45 m'i 



(7.41) 
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Substituting a value for QtXt (and including the colour degrees of freedom), this becomes 

,2 Q^7 71^2 in ^4 1 ^2 ][/r2 



F^'l''^ = Keg cos Ow 



(4^)2 

neg uo 
2cos Ow (47r)2 



11 m. 

14 + 



15 M2, 



37 Ml 19 m' 



3 M| 



1 — - sin^ ^vi/ 



15 M2, 210 M^, 210 M^, 
8 ^ 1 7m| + 11M| 



20 



45 



(7.42) 



and using My/ = Mz cos 9w, e = gsinOw and sin^ = 0.23, along with a numerical 
value for Mz/rrit |25|, we obtain the approximation 

,2 



(An) 



4 n 

18 + 1.7^+0.28-^^ 



(7.43) 



8. Calculation of the decay width 

8.1 Decays involving spin-2 Kaluza-Klein particles 

For the decay to a spin-2 KK excitation, the matrix element is given by 

^{h) ^ £^P*(q)e-*{k)ef^(p)F^f^) [{kxQu - k ■ qr^.x) {k^q^ - k ■ qr^^p) + (A ^ p)] . (8.1) 
To obtain the square of the matrix element, we need the polarisation sum formulae 



pol 




Y.e'^*{k)e^'{k) 



pol 



pol 



' Ml 



(8.2) 
(8.3) 

(8.4) 



(we note that the formula for the gravitational polarisation sum |^ |2^ differs from the 
massless case presented in Q). This gives [^] 



{Ml-ml)'{7Ml + 3ml) 



36M| 



(8.5) 



We note that as ^ 0, we recover the expression obtained for the case of decay to a 
massless graviton ]^], save that using the massive gravitational polarisation spin sum causes 
a difference by a factor of 7/2 from using the massless gravitational polarisation spin sum. 
The decay width T^'^'") to a single spin-2 KK mode is given by 



p(/i,n) 



327r2M| 



dn. 



(8.6) 
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where 



p* = [Ml - ml) 

^ 2Mz ^ 



(8.7) 



This gives 



Substituting the expression for F^^^ obtained in equation ( 7.20| ), and recalhng that a 
e^/Air and k = VSttG, we obtain 



2 \ 5 
n_ 

I 

0.00088 + 0.27 I 



1 \ 3 



EE 



(j + 2)(j + 3)(j+4) M| 



1 



i=0 j=0 



{i + 2){i + 3){i + 4)(j + 2)(j + 3)(j + 4) V M 



+ 



2 \ *+? 



To obtain the total decay width involving spin-2 graviton excitations, we must sum over 
the excitation levels. Following Q, we make a continuum approximation for the density of 
states at a given mass level. Given a common compactification radius R/2tt, the density 
of states over which we should integrate is 



n-2 



(47r)"/2r(n/2) 



(8.10) 



that is, the derivative with respect to m| of the volume of a n-dimensional hypersphere of 
radius r = That is, the total decay width is approximated by 



■ tot 



Ml 



.11 



We change the integration variable to mi/M^ and expand the resultant beta functions in 
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terms of gamma functions, obtaining 



r 



(h) 

tot 



727r2(47r)"/2 
0.00088 



7-5! 3- ^-5! 



r(f + 6) r(f + 7) 



+ 0-27 E 
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7-(j + 5)! , 3-f .(j + S)! 
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(j + 2)(i + 3)(i + 4) Vr(§+i + 6) r(H+j + 7) 



+ 



+ 21 < 



EE 

j=0 j=0 



{i + 2)(i + 3)(i + 4)(i + 2)0- + 3){j + 4) 



7-(i+j + 5)! , 3-^-(i + j + 5)! 



+ 



r(f + z + i + 6) r(f + i + j + 7) 



(8.12) 



We note that the width tends to zero as n — > oo, showing that as the number of dimensions 
increases for large n, the phase space increase dominates the increase in the number of KK 
states to reduce the decay width. 



8.2 Decays involving spin-0 Kaluza-Klein particles 

For the decay to a spin-0 KK excitation, the matrix element is given by 



Using equations (|8^ and (|8!^) , we obtain |]24| 



2\2 



(M|-m|) 



The decay width r(<^'") to a single spin-0 KK mode is given by 



■p((/),n) 



327r2 M| 



dn, 



where, as before. 



(8.13) 
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This gives 
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Substituting the expression for F^*^) obtained in equation ( [7.43D , again using that a = e^/47r 
and K = VSttG, and recalhng that uj"^ = 2/(3(n + 2)), we obtain 



a 



2887r2(n + 2)M| 

m?. mi nA 
330 + 63^ + 13^ + 0.97^ + 0.078^ 

M| Mi Ml Ml 



.18) 



To obtain the total decay width involving spin-0 graviton excitations, we again ap- 
proximate by an integral the sum over the excitation levels. Recalling that we have (n — 1) 
spin-0 particles at each mass level, the total decay width is approximated by 



r 



{4>) 

tot 



Ml 



.19) 



As in the spin-2 case, we change the integration variable to m^/M^ and expand the resul- 
tant beta functions in terms of gamma functions, obtaining 



,(^) _ g^GMf+'^fi" (n - 1) 
tot - 487r2(47r)"/2 (n + 2) 
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X 



r(§ + 4) 

0.078 • ^ 



63- f 

r(f + 5) 



13. |. (1 + 1) o.97-§-(f + l).(f + 2) 



r(f +6) 
(f + l)-(f + 2)-(f +3) 

r(f + 8) 



r(f + 7) 



.20) 



As in the spin-2 case, the decay width tends to zero as n ^ oo. We note that the width 
is zero for n = 1, which we should expect as there are no spin-0 excited states in the case 
of one extra dimension. We note also that whilst the width should be zero for n = 0, the 
expression above does not give zero for n = 0. However, the extra-dimensional spin sum 
identities of equation (f]B) are only valid for n > 0. 



8.3 Overall decay width to photon and Kaluza-KIein graviton excitation 



We may now use equations ( 8.12| ) and ( 8.20 ) to write down an approximation for the full 



decay width at one loop of the Z boson to a photon and a Kaluza-Klein excitation in the 
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ADD scenario: 



7-5! 3-f-5! 



0.00088 — + + 



r(f + 6) r(§ + 7), 
1 

(j + 2)(j + 3)(i + 4) vr(§+j + 6) ' r(f+j + 7) 



+ 27 <!V ^ f JMj±5)L + ^lUll^ M + 



1 + 2)(i + 3)(. + 4)(j + 2)(i + 3)(i + 4) 



7-(i + i + 5)! ^ 3-f -(i + j + S)! - ^ 



r(§ + i + j + 6) r(§+i + i + 7) 
3 (n-l) f 330 _ 63-^ 13 • H . (" + i) 



2(n + 2) lr(f + 4) r(f + 5) r(f + 6) 

0.97 • § • (§ + 1) • (f + 2) 0.078 • f • (§ + 1) • (f + 2) • (f + 3) 

r(H + 7) r(§ + 8) 

(8.21) 

Two comments about this expression are inherited from the constituent expressions of 
equations ( g.l2 ) and (|8.20| ). Firstly, the overah expression does not reduce to the case of 



a single, massless graviton for n = 0, because the extra-dimensional spin sum formula in 
equation ( |1.5D is only valid for n > 1. Secondly, the overall width tends to zero as n — > 00, 
so that this process becomes less distinguishable from a Standard Model background as 
the number of extra dimensions increases. 

Tables |^ and |^ give numerical values^ for the contributions inside the square brackets in 
equation ( ^.2lD . It is notable that the contributions from the decays to spin-0 excitations 
are dominant. 

In order to estimate the effects of the analytic approximations used for the integral 
J{X,Y,Z) in equation (7^), and for the integrals performed to obtain equations ( 7.2g| ) 



and (|7.40| ) , we have evaluated contributions to the decay widths by integrating numerically 
over the Feynman parameters and the KK mass-squared m^, using Mathematica. We find a 
deviation of approximately 5% between the analytically approximated and the numerically 
approximated answers for the decays to spin-2 excitations, and a deviation of less than 1% 



^Aside from the analytically approximated expressions for the <j) widths, which are calculated by hand, 
these are evaluated using Mathematica . The infinite sums are evaluated by truncation once the quoted 
level of significance has been achieved. The numerical approximations to the integrals are evaluated using 
the NIntegrate function with default options, and (for comparison) using a number of other numerical 
integration routines |2g] with their default options. 
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n 


h an. approx. 


h numerical approximation 


full 


nif 


W,t approx. 


2 


0.16 


0.16 


0.16 


0.16 


3 


0.056 


0.058 


0.057 


0.056 


4 


0.019 


0.020 


0.020 


0.019 


5 


0.0066 


0.0070 


0.0069 


0.0066 


6 


0.0020 


0.0023 


0.0023 


0.0022 



Table 1: Numerical values of the h contributions inside the square brackets of equa- 
tion ( |8.21 ) for some different values of n ('an. approx.' column), and the equivalent values 
when the contributions are calculated by numerical integration. The 'full' column shows 
the results of integrating numerically the full expression for the decay width (without mak- 
ing the approximations made in the analytic case) , and is the value with which the analytic 
expression is compared. The extra columns show the effect of performing numerical inte- 
gration under the extra approximations of neglecting the masses of the fermions other than 
the top ('?7i/ — > 0' column) and of assuming in addition to this that the W and top masses 
dominate the denominators of their respective integrands ('W,t approx.' column). The pur- 
pose of including such columns is to give an indication of how much these approximations 
change the final answer. 



n 


an. approx. 


(j) numerical approximation 


full 


mj - 





W,t approx. 


Num. change 


2 


5.3 


5.3 


5.3 




4.4 


5.4 


3 


3.9 


4.0 


3.9 




3.3 


4.0 


4 


2.2 


2.2 


2.2 




1.8 


2.2 


5 


1.0 


1.0 


1.0 




0.85 


1.1 


6 


0.46 


0.46 


0.46 




0.37 


0.47 



Table 2: Numerical values of the (j) contributions inside the square brackets of equa- 
tion ( |8.21 ) for some different values of n ('an. approx.' column), and the equivalent values 
when the contributions are calculated by numerical integration. Again the 'full' column 
gives the values with which those calculated analytically are compared. The extra columns 
show the effects of neglecting the masses of the fermions other than the top ('mj — > 0' 
column), of assuming in addition to this that the W and top masses dominate the denom- 
inators of their respective integrands ('W,t approx.' column), and of making this assump- 
tion but changing the numerators to eliminate M^r as done in the analytic calculation 
('Num. change' column; as noted in the text this is equivalent to considering higher order 
corrections in the series expansions for the fractions in the integrand). 



for the decays to spin-0 excitations^. We have also calculated numerical approximations 
subject to the assumptions of zero mass for fermions other than the top, and (additionally 

*The quoted deviations compare the values obtained from analytic approximation with the values ob- 
tained from evaluating the full integral numerically using Mathematica's NIntegrate routine. 
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to the previous assumption) of the W and top masses dominating denominators in their 
respective integrands. This indicates that the zero fermion mass assumptions are a very 
good approximation, but that assuming that the W and top masses are much greater 
than half the Z mass yields a poorer, although still acceptable, approximation. It is also 
notable that the numerical approximations of the integrals for the decay into (p suggest 
that the accuracy of the analytic calculation of the width is much improved by rewriting 
the numerators of the approximated fractions as we have. 



9. Conclusions 

We have evaluated amplitudes for Z-photon-Kaluza Klein graviton/gravi-scalar interaction, 
relevant to extra dimensional models in which the Standard Model fields are confined to a 
4-brane. 

In addition, we have evaluated in the ADD scenario a reasonable approximation to 
lowest order for the decay width of a Z boson to a photon and any Kaluza-Klein excitation 
of the graviton. This width gives an extra contribution to the channel of Z boson decay 
to photon plus missing energy compared with the Standard Model. The channel can 
provide stronger bounds on the compactification radius R/2'k for small numbers n of extra 
dimensions. 

The consideration of the decay channel involving spin-0 KK excitations (gravi-scalars) 
has proved significant: for the decays considered, these channels provide the larger contri- 
bution to the overall width. The consideration of processes involving the gravi-scalar has 
not been common in the literature. 

We expect our signal to be most significant when the centre of mass energy is equal 
to Mz-, as, for example, was the case at LEP. We therefore intend to calculate bounds 



upon the compactification radius from LEP data ||12|], given that the data appear to be 
in accordance with Standard Model predictions. Our process would be additional to the 
tree-level radiation of a spin-2 KK graviton tower from a photon. 

The amplitudes obtained for Z ^ hj and Z ^ (j)^ could also be used in the calculation 
of decays of a RSI KK mode into a Z boson and a photon. In RSI, the KK mode coupling 
is enhanced by a warp factor, making resonant production at colliders and interesting 
signature, worthy of study. Although h — > jZ and (j) — > 7Z are loop-induced processes and 
are therefore suppressed with respect to other, tree-level decays of a RSI KK mode, the 
channel would be useful for checking the couplings k of the excitations, a la reference |29|. 
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A. Dimensional regularization integral identities 



We write D = A — e, and in each case take the hmit as e ^ 0. 7 is the Euler-Mascheroni 
constant. The integrals on the left are in Minkowski space. The methodology used may be 



found, for example, in Appendix A. 4 of reference [30| 
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